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ABSTRACT 

The order dependent mapping method, its convergence has recently been proven 
for the energy eigenvalue of the anharmonic oscillator, is applied to re-sum the 
standard perturbation series for Stark effect of the hydrogen atom. We perform 
a numerical experiment up to the fiftieth order of the perturbation expansion. A 
simple mapping suggested by the analytic structure and the strong field behavior 
gives an excellent agreement with the exact value for an intermediate range of the 
electric field, 0.03 < E < 0.25. The imaginary part of the energy (the decay width) 
as well as the real part of the energy is reproduced from the standard perturbation 
series. 
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Quite often the perturbation series in quantum theory is not only divergent but 
even non-Borel summable [1]. Such a non-Borel summability is sometimes known 
to have its physical origin. For quantum mechanics with an unstable potential or 
degenerated potential minima the tunneling (barrier penetration) effect gives the 
Borel singularity [2]. For a uniform electric field in QED, it can explicitly be shown 
[3] that the instability of the vacuum state causes the non-Borel summability. One 
may even conjecture [3] that the Borel singularity in QED and QCD is an indication 
of the instability of the perturbative vacuum. 

The above observations in turn suggest a possibility that for certain case the 
standard perturbation series has enough information on the quantum tunneling (or 
the vacuum instability) that is usually regarded as "non-perturbative" effect. 

Recently Kleinert [4] has made an interesting observation that the decay width 
(imaginary part of the energy eigenvalue) for an unstable potential (anharmonic 
oscillator (AHO) with the negative coupling g < 0) can precisely be reproduced 
only by using the perturbation series of the energy eigenvalue of AHO, hence 
explicitly implemented the above possibility. His method is based on the so-called 
delta expansion [5] (variational perturbation method) that is a special case of the 
order dependent mapping method (ODM) [6,7]. The convergence of the method for 
AHO has later been rigorously proven [8] for \g\ > go (go is the radius of convergence 
of the strong coupling expansion go ~ 0.1). The proof itself [8] is applicable for 
other quantum systems if one has enough information on the analyticity and the 
strong coupling behavior of the interested quantity. 

In this article, we take a simple quantum mechanical system for which the 
perturbation expansion exhibits the Borel singularity on the positive real axes and 
apply ODM to re-sum the standard perturbative expansion. The physical origin 
of the non Borel summability in this system is the instability of the "vacuum." 

We consider the ground state energy of the hydrogen atom in a uniform electric 
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field (Stark effect), 

\-^—V 2 - - - eEz) = 5*. (1) 
y 2m e r J 

This system is unstable due to quantum tunneling (or spontaneous ionization) 
and thus the energy eigenvalue of the quasi-ground state has the imaginary part 

(=-r/2)*. 

On the other hand, the standard weak field expansion of the "ground state" 
energy reads ^ 

oo 

£(E)~J2a n E 2n , (2) 

n=0 

where the first several coefficients are 

a = -1/2, ai = -9/4, a 2 = -3 555/64 

a 3 = -2 512 779/512, a 4 = -13 012 777 803/16 384, (3) 
o 5 = -25 497 693 122 265/131 072. 

(An economical way to compute the higher order coefficients a n is summarized in 
Appendix.) The simple (truncated) sum of the series is of course real for E real. 
Note all the coefficients have the same sign. In what follows we will apply ODM 
[6,7] to this standard weak field expansion and see how the exact energy eigenvalue 
(that is complex) emerges. 

As was emphasized in [8] , knowledge on the analytical property and the strong 
coupling (field) behavior of the physical quantity are crucial for applying ODM. 

★ In what follows, the energy £ and the field strength E are measured respectively by, 2(Ry) = 
m e e 4 /h 2 = 27.21(eV) and l(a.u.) = m 2 e e 5 /h 4 = 5.142 x 10 9 (V/cm). Under a realistic 
condition in laboratory, the maximum strength of the electric field is E ~ 2 x 10~ 4 (a.u.), 
and practically is in the validity region of the standard perturbation theory and WKB 
approximation. The aim of this article is thus somewhat academic and to investigate how 
one can extract information from a badly divergent (and non-Borel summable) perturbation 
scries. 

f Since the perturbation hamiltonian has an odd parity, only the even powers of the electric 
field survives in the expansion for the ground state energy. 
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For the quasi-ground state of (1): i) The energy eigenvalue S(E) is analytic on the 
cut E plane at least in the vicinity of the origin and the cut is along positive real 
axis of E, as well as the negative axis, due to the quantum tunneling [9]. S(E) is 
analytic on cut — E 2 plane (at least in the vicinity of the origin), only the negative 
real axis being the cut. ii) The discontinuity along the cut is given by the WKB 
formula (for E > 0) 

lm£(E) = -Jy-Le^-^E-'e-V^il + 0(E)) ~ -2£~ V 2 /( 3S ). (4) 
1 + v2 

The above two properties in turn give the large order behavior of the perturbation 
coefficients in (2) [10] 

4 /3\ 2n+1 

an - -- (~ 2 ) (2n)!. (5) 

Corresponding to the instability of the ground state, the series is non-Borel summable. 
hi) The strong field behavior [10] is given by 

arg£ ~ -7r/3 + 0((lnE)- 1 ), 

o4/3 

\S\ ~ 2- 5 / 3 £ 2 / 3 (ln£) 2 / 3 + — ^(lnln^^nS)- 1 / 3 + 0(£ 2/3 (ln E)' 1 ^). 

o 

(6) 

This is an asymptotic expansion in contrast to the convergent strong coupling 
expansion in AHO case. 

In ODM [6,7], one first map the coupling constant E to another variable A by 
some function F(\) 

E = pF(X), F(A) = A + 0(A 2 ), (7) 
and factorize the original function with some prefactor function /(A) 

S(E) = /(A)*(A), (8) 

and expand ^/(A) with respect to A. The parameter of the mapping p is determined 
order by order by some condition. The choice of F(X) and /(A) is crucial for the 
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convergence property of the method, and is constrained by the analytic property 
and the strong field behavior of S(E) [6,7,8]. 

In the case of AHO [6,8], the (known) best choice of the mapping is equivalent 
to the delta expansion [5] which starts from a decomposition of the hamiltonian 



and the expansion with respect to 5 (and finally 5=1). The variational parameter 
Q (which corresponds to p in ODM) is determined order by order of the expansion 
by some condition. It can be shown [8] that the sequence obtained in this way 
converges to the exact energy eigenvalue for a wide range of complex g. 

As our first choice of the mapping in the present problem we may therefore 
start from a decomposed hamiltonian: 



where in the spirit of the delta expansion, we have introduced a variational pa- 
rameter Z. The expansion is done with respect to 8 (and finally 5 = 1) and Z is 
determined order by order by some condition. As for AHO [6], it is easy to see the 
expansion is equivalent to a mapping 



and a factorization 






(10) 



£(E) 



1 



*(A). 




(1-A)2 



The relation of parameters in both methods is given by (setting 5=1) 



The expansion of ^(A) up to iVth order reads 



^iv(A) = ]>> n (p)A 



n 



(14) 



n=0 



where 



bo(p) = ao, h(p) = -2ao, b 2 (p) = a + aip 2 , 



[n/2] 

bn(p) = £ 
m=l 




(15) 



As the condition to fix p order by order, we use [6,7], 



ttjv(A) - ttjv-i(A) = 0, 



(16) 



and pick up the real positive solution for p [6,8]. We may use the other condition 
such as [5] <9((1 — \)~ 2 ^ n) / 'dp = as well. In the case of AHO, the convergence 
property is known to be insensitive [8] to such a choice of the condition. 



At first sight, the choice in (11) and (12) seems to be well suited to the strong 
field behavior of S (6): For \E\ large, A ~ 1 and E ~ (1 — A)" 3 from (11), so 



logarithm in (6) cannot be reproduced in the present simple mapping. 

We have done a numerical experiment based on the mapping (11) and (12) 
to iV = 50. For an actual calculation of the perturbation coefficients a n , see 
Appendix. For a wide range of E, we observed a violating oscillating behavior 
such that ReS(E) > for E = 0.25 from the first several orders. The mapping 
(11) and (12) completely fails. 

The failure of the mapping (11) and (12) that is equivalent to a naive delta 
expansion like method (10), should clearly be related to the incompatibility of 
the mapping (11) and the analytic structure of £{E). As was shown for AHO [8] 
the mapping should be chosen to be compatible with the analytic structure. In 
this sense, the surprising success of the delta expansion for AHO was somewhat 



the prefactor in 



(12), (1 — A) 2 ~ E 2 ^ reproduces the correct power in (6). The 
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accidental. In ODM on the other hand we have a wide freedom to choose the form 
of the mapping and the decomposition*. In fact according to i) above the analogy 
with AHO case, for which the cut is only along the negative real axis, is not hold. 

This observation suggests a use of — E 2 instead of E as the fundamental variable 
for which the cut exists only along the negative real axis (at least in the vicinity 
of the origin). As the second choice therefor we take 

-E?=p—±—, (17) 

^(1-A) Q ' v ; 

and 

£ = 7 (18) 

(1-A)«/ 3 V ; K J 

The form of the prefactor in (18) has been chosen to reproduce the E 2 ^ behavior 
in (6). Here we again abandon to reproduce the correct logarithmic behavior in 
(6) and hence we expect the present mapping fails for a large E anyway [8]. The 
coefficients b n (p) in (14) now read 

KM = £ ( _ ir r( n + (^-i)m- a /3) 

[n — m)\Y{am — a/6) 



m=0 



For the value of a, we do not have any criteria at present (see below) and the 
choice is a sort of guesswork. We examined a = 3/2, 2, 5/2, 3, 7/2, 4 and 9/2. We 
have again done numerical experiments up to iV = 50. The p is determined again 
by the condition (16) and the real solution for p is taken. 

For a = 3/2, we observed a tendency of a slow convergence to a wrong answer, 
for example for E = 0.25, Re£ — > —0.63 and Im£ — > 0.1. For a = 2, the sequence 



* Of course it is quite possible that our choice of the variational parameter in (10) is not good 
enough to simulate the effect of the electric field. One may introduce more complicated form 
of the variational hamiltonian. Note however that for such a general choice, a connection 
with the standard perturbation series is lost and hence the advantage of the delta expansion 
is too. 
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again converges but to a wrong answer, namely Re£ — > —0.54 and Im£ — > 0.05. 
a = 4 and a = 9/2 cases do not work (slowly divergent). Among the above 
values, a = 5/2, a = 3, a = 7/2 seems work and a = 3 seems to give the fastest 
convergence. Therefore we will only report the result for a = 3 in the following. 

In Fig. 1, we have plotted the value of p determined in this way. A numerical 
fitting gives a scaling on the order N: 

p ~ 0.2851iy~ L686 . (20) 

In Figs. 2a and 2b, the relative error to the exact value of the energy (obtained 
by a numerical integration of the Schrodinger equation [11]) is depicted for a fixed 
value of the field strength E = 0.25 (2a for the real part and 2b for the imaginary 
part). A tendency of the fast convergence to N = 50 is observed. 

Observing this convergent behavior, in Figs. 3a and 3b, we plotted ReS(E) 
and Im£ (E) by ODM and the exact values [11] (full squares) in the intermediate 
range of the field strength. We see a surprisingly good agreement for the imaginary 
part as well as the real part. 

As was mentioned above, for a strong field regime, we expect the method fails 
because the present mapping can only reproduce the correct power behavior in 
(6) but not the logarithmic corrections. To see this we have plotted ReS(E) and 
Im£ (E) by ODM and the asymptotic behavior (6) (broken line) in Figs. 4a and 
4b. 

On the other hand, the experience in AHO case [8] suggest the method also fails 
to reproduce the imaginary part for the weak field region. In Table 1, the imaginary 
part obtained by WKB approximation (4) (which is a good approximation within 
this regime) and ODM are compared. It can be seen that ODM completely fails 
for E < 0.025*. 

★ For the real part of the energy, ODM gives an excellent result in this weak field regime. 
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In summary we have numerically observed that, at least for an intermediate 
range of the field strength, we can extract an accurate value of the imaginary part 
from the standard perturbation series by suitably choosing the mapping in ODM. 

In spite of this success, the convergence proof of [8] cannot directly be applied 
to the present case: First of all, a = 3 is a dangerous case in the view point of 
the proof [8] because the negative real of — E 2 is mapped to A with |A| < 1 which 
obstructs the direct application of the proof. One have to extend the "contour" [8] 
to the outside of |A| > 1 but this requires some knowledge on the analyticity on 
the higher Riemann sheet of — E 2 that is lacking at present. 

Secondly the prefactor (1 — A) -1 (for a = 3) in (18) is not sufficient to make 
^(A) finite as A — > 1, because of the logarithmic correction in (6). The boundedness 
of ^(A) at A = 1 was another necessity condition for the convergence proof in [8]. 

Finally the fact that the observed power in (20) is larger than —2 is rather 
strange. Repeating the analysis of [8], the error of the method caused by the 
divergence of the perturbation series (R^ part in [8]) is estimated as 



which diverges if the scaling (20) retaines for N — > oo and is in fact ~ 156 for 
N = 50 and E = 0.25 (The relative error ~ 26300%). This clearly contradicts with 
the observed fast convergence in Figs 2a and 2b. The convergence mechanism in 
[8] therefore cannot be applied in the present case. We need some new explanation 
such as a cancellation among R^ part and R^ part in [8]. The fact that ^(A) 
is not bounded at A = 1 also suggests this possibility. 

Obviously more theoretical study is needed to explain the empirical success 
of the mapping (17) with a — 3. Nevertheless we believe our observation gives 
another concrete support for the possibility that one can in practice reproduce the 
exact property (tunneling or decay in the present case) of the physical quantity 
only relying on the standard perturbation series. 
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It is possible to reproduce the correct logarithmic behavior in (6) by using more 
complicated mapping such as the one in [7] . Although we do not try in this article, 
it is crucial if one wants an accurate value for a large E. 

We thank Prof. T. Amano for useful information. The work of H.S. is sup- 
ported in part by Monbusho Grant-in-Aid Scientific Research No. 07740199 and 
No. 07304029. 



APPENDIX 

It is well known [12] that the Schrodinger equation (1) has a separated form 
in the parabolic coordinate x = y/Xpcosip, y = ^O^psimp and z = (A + a0/2. 
The problem is then reduced to one-dimensional and we may apply the recursion 
formula technique in [13]. Setting 

* = (Zv) m/2 e- { t +v)/2 F(Z)G(Tj), (Al) 

(m is the magnetic quantum number) where 

f = v /= 2£A, T) = v /= 2£>, (A.2) 



we have 



iF" + (m + 1 - OF' - Qm + \ - a - ^ F = 0, 
ZG" + (m + l-Z)G'- Qm + i-T+|e 2 )G = 0. 



(A3) 



Here a and r are the separation constants and the "coupling constant" g is defined 
by 

° +T= vki' 9= (-^w-- {AA) 

In what follows we consider only the case m — 0. 



10 



We define the "perturbative expansion" 

oo oo 

* = $>*<A F(0 = J2 F ^)9 J - (A.5) 

k=0 j=0 

The zeroth order solution is given by ao — 1/2 and -Fo(0 = 1. The equation (A. 3) 
has a polynomial solution order by order in g. Hence we further set 

m) = J2hi^ fo r j> 1 - ( A -6) 

Eq. (A. 3) then becomes a set of recursion relations: 

<rj = (A.7) 

and 

j_1 1 

(Z + l) 2 f h i +1 - lf h i + °*h-H - = 0, (A.8) 

k=i 

for j > 1 (/ ,o = 1 and / 0) , = for I > 0). 

We have numerically solved (A.7) and (A.8) by FORTRAN in double precision 
up to CTioo- We then have used a symbolic manipulation package MATHEMAT- 
ICA to invert (A. 4) (note r = YlT=o(~ ^) k(T k9 k ) f° r obtaining the perturbation 
coefficients a n in (2) up to 050- 

To examine the quality of our numerical value for a n , we plotted the ratio of 
a n to the asymptotic behavior (5) in Fig. 5. 
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Table Caption 

Table 1 Comperison of ODM with WKB formula (4) in the weak field region. 

Figure Captions 

Fig. 1 Scaling of p with respect to the order of ODM N . p is fixed by the condition 
(16). 

Fig. 2 The relative error of ODM (17) with a = 3 (2a for the real part, 2b for the 
imaginary part) with respect to the order N. 

Fig. 3 The real part (3a) and the imaginary part (3b) of the energy obtained by 
ODM in the intermediate region of the field strength. The full squares are 
the exact values [11]. 

Fig. 4 The real part (4a) and the imaginary part (4b) of the energy obtained by 
ODM in the strong region of the field strength. The asymptotic behavior (6) 
is depicted by broken lines. 

Fig. 5 The ratio of the perturbative coefficients a n and the large order behavior (5). 
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E 


Im£ (WKB) 




Im£ (ODM) 




0.01 


-2.44847 x 10- 


-27 


-1.38778 x 10- 


-16 


0.0125 


-1.20942 x 10- 


-21 


+4.16334 x 10" 


-16 


0.015 


-7.30789 x 10- 


-18 


+1.09635 x 10- 


-14 


0.0175 


-3.58320 x 10- 


-15 


+7.32747 x 10- 


-14 


0.02 


-3.66731 x 10- 


-13 


+3.27433 x 10- 


-13 


0.0225 


-1.32341 x 10- 


-11 


-1.04966 x 10- 


-11 


0.025 


-2.30534 x 10" 


-10 


-1.66882 x 10- 


-10 


0.0275 


-2.36688 x 10- 


-9 


-1.66460 x 10- 


-9 


0.03 


-1.63588 x 10- 


-8 


-1.11321 x 10- 


-8 



Table. 1 
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